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CONVERGENCE OF GENERALIZED MUSCL SCHEMES 


STANLEY OSHER 


Abstract. Semi-discrete generalizations of the second order exten- 
sion of Godunov's scheme, known as the MUSCL scheme, are constructed, 
starting with any three point "E" scheme. They are used to approxi- 
mate scalar conservation laws in one space dimension. For convex con- 
servation laws, each member of a wide class is proven to be a convergent 
approximation to the correct physical solution. Comparison with another 
class of high resolution convergent schemes is made. 


^1. ^Introductk^^ Recently, there has been an 

enormous amount of activity related to the construction and analysis of 
"high resolution" schemes approximating hyperbolic systems of conserva- 
tion laws. Some examples of the successful consequences of this 
activity can be found in the proceedings of the latest (sixth) AIAA Com- 
putational Fluid Dynamics Conferences [1], CIO], [18]. Extensive 
bibliographies can also be found in these papers. 

Our aim here is merely to construct and prove the convergence of a 

subclass of these schemes approximating scalar convex conservation laws. 
# 

This subclass is based on an idea of van Leer [17]. He christened his 
algorithms MUSCL (Monotonic Upstream Centered Schemes for Conservation 
Laws) schemes, and with mixed emotions, we shall use his acronym here. 


In future work with S. Chakravarthy , we shall extend this construc- 
tion to systems in multi-dimensions, using triangle-based algorithms. 
This work in progress will stress the computational aspects of the algo- 
rithms, especially as they relate to the Euler equations of compressible 
gas dynamics. In earlier work with the same author, [123, we con- 
structed, and proved convergence of, a class of high resolution schemes 
approximating scalar convex conservation laws. We also showed for cer- 
tain high resolution approximations to systems, that limit solutions 
satisfy an entropy inequality. 

We shall consider numerical aproximations to the initial value 
problem for a single conservation law in one space dimension 

( 1. 1) (a) w fc + f(w) x = 0, t > 0, -1 < 1, 

with a periodic boundary condition: 

(b) w(x + 1 , t ) - w(x,t) 
and initia?. condition: 

(c) w(x, 0) = v/qv'x). ' 

It is well-known that solutions of ( 1 . 1 ) may develop discontinui- 
ties in finite time, even when the initial data are smooth. Because of 
this, we seek a weak solution of (1.1), i.e. a bounded measurable func- 
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tion w, such that, for all <S 4 C Q (R x R* ) , 

(1.2) (a) tS (w<t. + f(w)d) dx dt = 0 

RxR + 1 x 

(b) lim ||w(x,t) -w n (x)M * = 0. 

t-»0 L 1 

Solutions of (1.2) are not necessarily unique. For physical rea- 
sons, the limit solution of the viscous equation, as viscosity tends to 
zero, is sought. In the scalar case, this solution must satisfy, for 
4 Cq°^R x R + ), i >_ 0, and all real constants c: 

(1.3) (a) -» XX (|w - c !d fc + sgn(w - c)(f(w) - f(c))cf x ) dx dt £ 0 
This is equivalent to the statement: 

(1.3) (b ) |w -c| + ^ ( ( f ( w ) - f(c)) sgn(w - c)) < 0, 

in the sense of distributions. 

Such solutions are called entropy solutions. Kruzkov has shown in 
[93, that two entropy solutions satisfy 


(I. 1 *) I iw(x,t 1 ) - u(x,t ^ ) | 

for all t 1 t Q . Hence, (1.3) 
to the scalar version of (1.2). 


I t < ! lw(x,t Q ) - u(x,t Q ) | | 1 
L L 

guarantees the uniqueness of solutions 
Existence wa3 also obtained there. 
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Equations (1.3) may be viewed as the statement 

(1.5) V(w) ^ F (w) < 0 

for each of the convex "entropy” functions of solutions w to (1.1): 

V(w) = |w - c! 

and their associated flux functions 

F(w) = (f(w) - f(c)) sgn(w - c). 

For scalar convex conservation laws whose solutions lie in the space 
BV, a single entropy inequality (1.5) for any strictly convex V(w) 
and its associated entropy flux F(w), satisfying 

(1.6) F»(w> = V'(w)f*(w) 

is enough to imply existence and uniqueness of solutions to (1.1). This 
follows from the results of DiPerna [33. 

Next we consider a semi-discrete, method of lines, approximation to 
(1.1). We break the interval (-1,1) into subintervals: 

Ij = ( X 1 ( j - V 2 ) & < X < (j + V 2 ) & } 

j - 0, ±1 ;tN, with (2N + 1) A = 2. 


M 


1 t 
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Let Xj = j & , be the center of each interval 1^, with end 
poi.ts Xj.1 / 2 .x j+1/2 . 

Define the step function for each t > 0, as 

U^(x,t) s Uj(t), 

for *41., 

J 

The initial data is discretized via the averaging operator T^ f 

T /( ^ g (x) = ^ w Q (s) ds = Uj(0) for x 4 Ij. 

J 

For any step function, we define the difference operators 

*Vj 5 ±(u jti - u j> 

D ± u j * k /A ± u j 

A method of lines, conservation form, discretization of (1.1), is a sys- 
tem of differential equations 


( 1 . 8 ) 


SF u j * D + h j-i/ 2 s J = °» ±1» •••» ±N 


U^(x,0) = T /N w Q (x), for x 4 Ij. 


Here, the numerical flux defined by: 
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0.9) h j-V 2 5 h ^ u j+k-l’ * * * * u j-k^' 

for k > 1, is a Lipsohitz continuous function of its arguments, satis 
fies the consistency condition: 

h(w,w,. . ,w) = f(w). 

It is well known that bounded a.e, limits, as A 0, of approxi 
mate solutions converge to weak solutions of 0.1), i.e. (1.2) (a) is 
satisfied. However, this does not also imply that the limit solutions 
will satisfy any of the entropy conditions (1.5), let alone the general 
condition 0.3). Some restrictions on h are required. 

A simple class of flux functions h, for which (1.8) converges, 
for all f, to the unique entropy solution in L oo (L 1 (R);[0,T]), as 
> 0, for any T > 0, is the class of n E" schemes introduced in 
[11]. Such schemes satisfy the following: 

A consistent scheme whose numerical flux satisfies 

OOO) sgn(Uj - u j_i>E h j, -v 2 - f(u > ] i 0 

for all u between u j_i and u j is said to be an E scheme. This i.s 
at present, the most general class of schemes known to converge in the 
nonconvex case. 

It is clear that this class includes the widely known class of 
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three point monotone schemes, i.e. those for which 


h j-V 2 5 h(u j ,u j-i ) 


with h. i, noninoreasing in its first argument, nondecreasing in its 


second . 


We denote partial derivatives of a numerical flux via 


bu j+y h(u j + k ,,, * ,u j-k-f1 ) = h r 


Thus a three point scheme with a differentiable flux function is mono- 
tone iff 


h^ < 0 £ h Q . 


Any numerical flux can be written 


(1 - ,1) V/2 ■ V 2 (f< V * f<U j + 1 >) - teW AW 

where Q. i. can be viewed as the viscosity of the scheme [15]. 


One particular three point scheme is due to Godunov, [5], and has si 
special significance in this theory. The flux for Godunov’s scalar 
scheme can be defined by 


(1.12) IT -|. = h p (u .,u . .) = min f(u), if u, , < u, 
j-'/ 2 G j j-1 „ j-1 ~ j 


U j-1^ U j 
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max f(u), if u i 1 > u 


One can thus, [11j, characterize E schemes as precisely those for 


which: 


(1.13) (c) 

h .1-’/2 - h J~V if “j < U J-1 

(b) 

h j-V 2 ^ h j-V lf VVr 


or, as Tadmor [15] pointed out, those which have at least as much 
viscosity as Godunov's scheme 


<1 - 14) 9 j-V 2 * VV 

It follows from [11], Lemma (2.1), that these approximations are, 
at most, first order accurate. 

Together with an entropy inequality, a key estimate involved in 
many convergence proofs, is a bound on the variation. Fo," any fixed 
t ^ 0, the x variation of U^(x,t) is 

B(U^) = Z IA + Uj(t)|. 

J 

If we can write 


( 1 . 15) (a) 



A + u j ♦ D ._ 1/2 £ 
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lDJ V/2 i 0 

(0) V / 2 i 0 

then it is easy to show, [12], using an argument of [19], that, for 

S > t 2 > o 


B(U^C,t 1 )) < B(U^(*,t 2 )). 

Harten in [7], pointed out for explicit methods, that this decompo- 
sition could be obtained for schemes which are higher order accurate. In 
our present method of lines context, it involves a five point, con- 
sistent, approximation: 


(1.17) (a) 
with 


H , 

« = ' £ A * h i-'/ 2 = a S*V 2 & + U J " D j-V 2 A _“j 


c i*y 2 = c< V' u j*i’ u r"H> i 0 
<0) d J- 1 / 2 = D(u j+i ,u j ,u j-r u j- 2 ) i °> 

both Llpschitz continuous functions of their arguments. (See also van 
Leer [16].) 

In addition to (1.16) we have a maximum principle for (1.17), 

[ 12 ],: 


k 


f: 


10 


V 


t 
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(1.18) min u k (0) < Uj(t) < max u }< (0) 

k k 

for eaoh j and all t 0. 

Morover, in [12] we also showed a limit on th* possible accuracy of 
approximations of the type (1.17). Any scheme of the type (1.17) is at 
most first order accurate at nonsonic critical points of u. (A sonic 
point u is one such that f ' (u) = 0). Thus, although schemes of thi3 
type can be made to be as high as third order accurate, Lipschitz con- 
tinuity implies a local degeneracy to first order accuracy at smooth 
maxima and minima. This local degeneracy, together with some results on 
initial boundary value problems in [6], indicate strongly that overall 
second order accuracy is the best possible. 

Following Harten [73, we call algorithms of the type (1.17), total 
variation diminishing, or TVD schemes. 

It is known, that although (1.17) generates a compact family of 
solutions, U^(x,t), in L^U 1 (R) , [0,T]) , (if w Q L 1 f| L 00 fl BV) ), 
limit solutions may not satisfy any of the entropy conditions (1.15), 

and hence may not even be unique (1.12) See e.g. [4], [11], and [153. 

> 

One way of avoiding this difficulty is to use an inequality obtained in 
[113., 

Let V(w) be any convex function. We showed in [11], Section III, 
that, for any solution of any scheme (1.8)? 


11 — 
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(1.19) A(^ V( Uj > + D^CUj)) = X^ +1 dw V»(w)Ch j+1/2 - f(w)3, 


where the approximate entropy flux is defined through 


(1.20) F A ( Uj ) = F( Uj ) + VMUjUh^y - f( Uj )3. 

Thus, a sufficient condition that any limit solution satisfy (1.5) 
for a fixed convex V, is that 


( 1 . 21 ) 


>1 

'j 


dw V"(w)[h 


J+V2 


- f(w>3 


In order that the above inequality be valid for all convex V, it 

is necessary arid sufficient that h. correspond to an E scheme, 

J +/2 

which implies again that the approximation be at most first order accu- 
rate. Thus, for the schemes to be constructed in the following section 
we shall only obtain our entropy inequality (1.21) for a single V(w), 
say V(w) = V 2 W 2 . The following Theorem summarizes the technical 
hypotheses needed for convergence. 


THEOREM 1.1. The sequence of approximate solutions converges a.e. 
as 0, t£ the unique solution of the scalar , convex conservation law 

(I*!) Provided that the initial data ijs in BV and that inequalities 
(l.JX) and (1.21) are valid for £ single convex V(w). 


The proof is analogous to that of Theorem (4.1) in [123. 
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Xl^ Co n 3 1 r u o t to n of T V D MU 3 CL So h eme 3 . The MUSCL scheme as dis- 
cussed e.g., in [23 and [83 is a second order accurate extension of 
Godunov's method [53 that is based on ideas expressed by van Leer in 
[163. 


We assume that u^(t) is known. The first step in MUSCL is the 
reconstruction of a piecewise linear description of the solution. In 
the interval lj the result of this operation is 

(2.1) z (x i t ) = Uj(t) + (x - Xj)sj(t), for (x - x^) < A /2. 

Here Sj(t) are slopes that satisfy 

Sj = w x (Xj,t) + 0( A ) 

subject to some monotonicity constraints discussed below. 

Godunov's scheme, in our semi-discrete context, is the following. 

n 

The numerical flux h (u^.|,Uj) is computed by solving the Riemann 
problem, i.e. the initial value problem (1.1) with initial data 

w - Uj(t) , for x < x J+1/2 

* 

w r u j+1 (t), for x>* j+ 1 /2 . 

The resulting unique entropy condition satisfying solution is a 
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function of the type 


* 
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= w (£). 

A closed form for w was recently obtained in [11]. 

Then Godunov’s flux is defined through 

(2.3) h° +1/ * f (w(0) ) 

which is the same as (1.12). 

A fully discrete, explicit in time, Godunov scheme has the same 
numerical flux, after we impose a CFL restriction which prevents the 
interaction of solutions of adjacent Riemann problems The literal 
second order accurate extension to this fully discrete approach would be 
to compute the solution to the initial value problem (1.1), with initial 
data 



w - Uj(t) + (x - Xj)Sj(t), for x < *j + l /p 

- u j+1 (t) + (X - Xj)3 j+1 (t), for x > x J+1/2 . 

Find f(W(x j + i^,s) ) for t + ) s ;> t, with \ > 0 satisfying a CFL 
restriction, then compute 

(2.4) 5 ^ f(w<x jtV2> s» ds. 
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This is the MUSCL numerical flux. 


a 


! 
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Unfortunately, obtaining the exact solution to this nonlinear ini- 
tial value problem with piecewise linear initial data is a non-trivial 
business, even in the scalar case. However, at s = t + , x = x..j, , the 

J+'c 

solution is the same as for the Riemann problem with initial data: 
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w r Uj (t ) + ^Sj(t), for x < Xj + 1/ 2 

' ! vzy Ui for 

which is easily calculated, e.g., from the formula in [11], Thus, the 
semi-discrete MU5CL extension of Godunov’s algorithm comes through 


(2 -5> -St 1 ■ &- hG < u J+ i s j+r u j * 2 s j>- 

We generalize this still further as follows. Let 
h(Uj + ^,Uj) = toj + 1/ 2 he sn arbitrary first order accurate numerical flux 
function. Then our generalized MUSCL algorithm is merely: 


( 2 . 6 ) 


& 


- A ^- h(u j+ i -f 3 j*r u j + f 3 j>- 


We have: 


LEMMA 2.1. At_ points w(x^,t) = u^(t) i£ £ neighborhood of which 

p 

h(u. „,u.) is C with Lipschitz continuous partial derivatives, and 

J+1 j • — 

, the algorithm (2.6) is at least second order 


A 3 i4.l As, 

= 1 ♦ o<A) - q-J 


j 


15 


1 


accurate for smooth functions w. 
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Proof. Let 


(2.7) 


“j+V 2 = V 2 ( %1 + U J>- 


We shall show: 


(2.8) A. (h(u j+1 »J> - f <'V/ 2 ,) = o'A 2 ) 

at these points. 

Using consistency of h(Uj +1 ,Uj), the left side of (2.8) can be 
written as 


A 


-M u 






1 .U, + 


^ u 


j 


J T 2 "J 


1 - 


& s 


J±l! 


A u 


Jj 


w 0 




dy 






A s . 
J 


^ u jJ 


- 1 


The result is immediate. 


Next we make an observation about the viscosity of this scheme, 
, as defined in (1.11). 


LEMMA 2.2, Lf h ij3 monotone , then the viscosity of algorithm (2.J5) 

s ui s. 

is*a decreasing function of ’ and •, d . 

s /^ Uj — 2^-07- 
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The proof comes directly from (1.11). 

a i 3 i+1 

Thus, as expected, if and — are both restricted to be 

/ \ U j f \ J i 

j -* j 

positive, the most viscous, hence the least accurate, case occurs when 
these values ars both zero. The scheme then degenerates to first order 
accuracy. We can compress or smear out the solution locally by increas- 
ing or decreasing these ratios. 

We wish this scheme to be TVD, which will give a natural restric- 
tions on these above ratios. 


LEMMA 2. 3. ^ h h a flux corresponding to an E scheme , then the 

scheme is TVD if 


(2.9) 


il 3 < s 


01 / -V u j ’ £v u 


j±i 


< 1 


for each j . 


Proof. We have 


(2.10) -fsj.Uj., + 


- th(u -A 


A 


A L,, ">1" 2 “ ,u j + T s j ) “ f(u j + T s j ) 


A 3 1 A 

+ f(Uj + “ h<Mj - — ^Au^ +^Sj)] 
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+ i Ch(U j “^ U j + ^ S j } - f(U j “ 2 ) 


^ S 1 ^ A 

f(Uj - — gT*) - h(Uj 2 ' u j-1 2 3 j-1 


^3 


)3 


Thus, in (1.15) we can write: 


J±1 


<2. 11 Ha) C j+1/2 KVl-t^“i + t 3 J> - f(u j + 7 3 j 


) 


f(u J +l 2 3 J> - htu j -^ S j' U j + ^ 3 j >] 


(b) D. 


>v 2 = S"ST5jr Ch(u j -“ ,u j + - 2 °j 


&3 


i, U, 4-^3, > - f(4. -£s.) 


f(U j -f S j ) - h(U j 'f P j' U j-1 + ^ 5 j-1 )3 

Both ^j + 1/ 2 » D j 1/ 2 are nonnegati.ve because of (1.10) and (2.9). 

The restrictions (2.9) can be relaxed somewhat if h is monotone; 
we have: 


LEMMA 2.M. If h is a flux corresponding . to a monotone scheme , then 
the scheme is TVD if 


(s. - - s.) & 

(2. 12) (a) 1 > - -^~2~7\ — u ~ — near points where h^ 0 
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i s it! J — - near points where h 0 i 0 

2/\_Uj u 


for each j . 


Proof. We have 


(2.13) h(Uj - £ Sj.Uj., + f "j.,) 


As 


itl 


= A th(i vi 2 —'“J 


u, + 


as, , a^ a 1 ),, 

J-) - h(u, - — 2 * u j + “"" Q~' ■* 


As. &s. a*. ^WL- 

+ A th(U j " “2^’“^ + _ 2 J ‘ ) * hlU 3 " * 2“ 


A=. 


As,.. 


Thus, in (1.15). we can write: 


(2. 14) (a) C 


j+V? 


aroVVi + ^ )[/ A u j -$*j*i ' V 1 


(b) 


and 


Vv 2 - a<AjJT h o < u 

Both C >1/2 and D >1/s 

(2.12). 


j •■^r L ' s j )t/ -v u j + ^ (s j - 5 j-i )] 

are nonnegative because of monotonicity 


Remark 2.5. These last two Lemmas can be made "local," i.e. in 
regions of monotonicity we may relax the restrictions from (2.9) to 
( 2 . 12 ). 
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III. ,Co^n 5 en^^ In order to P rov<! convergence 

for convex f , we need only verify the discrete entropy inequality 
(1.21) for the schemes constructed in the previous section (Theorem 
(1.1)). We shall do this for the entropy 

V(w) = V 2 w 2 

and for schemes satisfying the TVD hypothesis of Lemma (2.3)» in addi- 
tion to some other restrictions near sonic shock points. 


Suppose Uj < Uj +1 , i.e. we are at a rarefaction. We then have, 


letting u j + i/ 2 = V 2 (Uj + u^ +1 ) as a^cve 


u i+1 As, . As, 

(3.1) 4^ dw[h(u j + i " 2 ,u j + 2 


,u, + - ■-*- ) - f (w) 3 


As 


j+1 


As, 


■ V (h( %i - u j + -r> - f( v/ 2 » 


u . 


♦ * U ; +1 dw [ f ( u j + v 2 > - «■» 


= ci] + [in. 


Now, using the hypotheses of Lemma (2.3)» we have: 


(3.2) 


[I] < 0 
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(3.3) (ID = da f“- 1+1/2 [^-(d - w)j f'(s) 


= da (u. - a) f'(u. i, ) 


- S ^ +1 dw / (s - w) f ** ( s ) ds 


’ ^ +1 dw / ^ /2 (s - w) f"(s) ds 


= - S J ' dw / 

U j 


- 1/2 Jy 1 ds f”(s> (minis - u^ +1 l,|s - u^ir < 0. 
J 


Suppose u^. > Uj +1 , i.e. we are at a shock. Let a j + 1/ 2 


be chosen 


so that 


4V, j+1 f'(w)(w - a j+1/2 > dw = 0 
J 


( 3 . 4 ) 


/J + 1 «f(»)dw 

ff j +, / 2 ' u = ^v f(u J > 

V f ' (w) ^ J 

j 
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where F(u) Is the entropy flux corresponding to our entropy, 

V(u) a Vgu . Of course we need the denominator in (3.4) to be nonvan- 
ishing. In fact we restrict ourselves to the situation where 


(3.5) 


We note 


U j+1 ~ a j+V 2 - V 


(3.6) 




'j+1 


(w - u j + i/ 2 ) f’(w) dw 


j+i 


(w " U j+V 2 ) - f, < u 


J+Vfc 


)) dw 


/\ fCUj) 


Thus, we have 


sup f»(u) 


(3.7) 


1 




' V/2 ” “j+tfe 1 < 6 ^ ,J j ) l» j+ v 2 5 


f<u.) 

where s. h = =■*— . 

j+V 2 /\ + Uj 


We further restrict the s^ so that 


(3.8) 


0 <. - A s j < 2 max(min((Uj - 3 j + 1/ 2 ) , ( 3 jJ / 2 " Uj)),0). 


In view of (3.7), this restriction does not affect the second order 
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aocuraoy of the soheme, except at sonic points. 


We now write 


(3.9)/ U , J+1 dw[h 


A S ■< + <| A, 3 i I 

f—t U 4 + O - ?<*)] 


j + 1 


J 2 


u j [ti 


U J+1 2 ’ U j + 2 




d « tf <»> - f( v/ 2 >] 


= [III] + CIV3. 

Now, the assumptions (3.5), (3.8) imply 

(3.10) till] < 0. 

An argument analagous to that of (3.3) suffices to show that 

im < o. 

Finally, suppose u^ +1 < u^ and (3.5) is violated. The inequali- 
ties of Lemma (2.1), and (3-8) imply that Sj and s^ +1 vanish. Then 
inequality (1.21) is immediate, since h is an E flux. 


We may summarize all this in 
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THEOREM 3.1. The sequence of approximate solutions satisfying 
(2.6) converges a. e. to the, unique solution of the scalar oonvex conser- 
vation law (J,. j.) provided that the initial data is in B.V. and that for 
e ach j 


0 < 7T~!” < 1 




i 


A ® 4 

o < Tr-e 1 - < i 




J 


and , in addition if u^ > u j + -|i 


“ A s j < 2 max <roin (< u j - a j-1/ 2 “ u j)>* °). 

VL ^Steady^Discret e^Shocks. We now check for the existence of 
discrete, steady, shock solutions to the convergent MUSCL schemes con- 
structed in the previous section, based on the Engquist-Osher flux, [4], 
Since the scheme satisfies a maximum principle (Lemma (2.1) of [12]), 
any profile must be monotone. 

Let f"(u) > 0, and for simplicity, we take f‘(0) = 0 = f(0). In 
this case, the monotone E-0 scheme becomes 


( 4. 1 )(a) h E0 (u j+1 ,u j ) = fju j+1 ) + f + (Uj) 

where 
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<b) 

Mu) 

= f(u) 

if 

u < 0 


Mu) 

5 0 

if 

u > 0 

Co) 

f + (u) 

s f(u> 

if 

u > 0 


V»> 

r 0 

if 

u < 0 


I R 

Let u , u be the left and right states for a physically correct, 
steady shook, i.e. 


f(u L ) = f(u R ), u L > 0 > u R . 


A steady discrete shock (u.) will satisfy 


(4.2)(a) lira u, = u 


j ~oo 


J 


(b) lim u. s u 


R 


j-iM-00 


j 


(q) h 


EO 


£P i+ i A 5 ii L 

U j+1 “2 9 U j * 2 s f ( u ) 


= f 


'J+1 



/ 


+ f. 



for all j. 
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We shall seek steady, discrete, shocks of the same general form 


obtained in Cl] for the E-0 scheme, o.f. also [12], section V. 


These are 


(4.3)<a) 


u j - u » j < -1 


Uj - u n , j > 2 


uL > u q > 0 > u 1 > u R . 


For the first order scheme, u Q can be viewed as a smooth function 
of u^ satisfying the above, and, in addition 


f(u Q ) + f( Ul ) = f(u L ) 


For the present scheme, we shall get a different one parameter family of 


intermediate states. 


It follows for all j ^ 0, that (4.3) implies (4.2)(c). 


For j = 0, we have the following equation: 


3i & I 

(4.4) 0 = f_ u, - + f + u Q 4- - f(u L ) = G(u 1 ,u 0 ) 


One special solution is, again, = u , u Q = 0, 


u 

For (Uq,u^) close to (0,u ) and satisfying (4.3)(c), we have, 


m 
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from the hypotheses of Theorem (3.1): 


(*».5) a, = 


(u - u,) 


1 ■- fl V 2 A 


>o 5 0 


Here 0 £ 3 1 a is assumed, in addition, to be a C * function of 

bai 


near = u , with 


r 2 - s 0. 


A straightforward calculation gives U 3 


C 4 . 6 )(a> ( u R ,o) s f»(u R ) 


i + -j£l < o 


(b) jrp (u R ,0) s 0 
0 


(c) ^<u R ,0) 


bu 




■*2 


a 1 I 

1 + *rj >0 


(d) <uR -°) ■ 0 


v 2p 

(e) — = f»(0) > 0 

»“o 


A simple application of the implicit function theorem gives us: 
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THEOREM 4.1. Given the hypotheses of Theorem ( 3 ,. j.) the mild 

technical •SB<!rnl!!& ’1/ 2 ' a - ^ SOSb 

discrete , shock solutions to (2. 6) of the type (4. 3.) , with u 1 a smooth 
function of u Q , and 0 < u Q , small enough . By sy mmetry , the same is 
true , with u Q a smooth function of u^ < 0, and — u ^ small e nough . 


V. Comparison With Other Convergent High Resolution Schemes. In 
[12] we constructed a class of convergent high resolution schemes, using 
flux limiters and second order accurate upwinding. Our first step was 
to use a first order E scheme flux, h(Uj^,Uj), to construct a second 
order accurate TVD scheme. The resulting scheme is 


i 

*| 

i 

i 


: i 

V 

*, 

t 

i 


(5.1) 


<5u 


j 1 

"St = “ /\^- H(u j+2, u j+r V u j-i ) * 


with 


(5.2) H(u j+2 , u j+r u jf u^) a h<u j+r u^ - V 2 * (R‘ +1 )(h(u j+1 , Uj ) - f(u.)) 

+ \ l KRj)( f ( u j +1 ) - Mu J+1 , u j)>* 


N 

n 


Here 


(5.3) 


f(u.) - h(u, , u. .) 

'n^VV 


h<u j + i ,IJ .i > - f(u i ) 

3 "^y Vi> - f( Vi> 


' £ 

4 

f 

f 

1 

'f 
* 9 . 
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The slope limiter i) (R) satisfies i)(1) = 1 and the following 
inequalities 


(5.4)<a> 


1 + V 2 


* (HT) 

— L - * < R j + i> 


R. 


> 0 


Cb) 


1 + V 2 


* (Rt> 


41 - 




> 0 


Various slope limiters have been developed. See Sweby [133 ( for a 
numerical and theoretical analysis of their properties. Perhaps the 
simplest example i3 the so called min-mod, [14]: 


(5.5) 


i)(R) = max (0,min(R f 1 ) ) . 


We have the following analogue of Lemma (2.2): 

LEMMA 5.1. If h is an E flux , then the viscosity of (5.1.) is. £ 
decreasing function of i?. 

Next we compare the schemes (2.6) and (5.1). If f'(w) i 0 and w 
is smooth at Xj, then it is easy to see that 


H(u j+2’ u j+1» V u j-1* = 


^j-1 ^1! 

u j + 1 " 2 ’ U j + T 


+ 0(QX Uj) 2 + (/^ Uj) 2 ) 


if 
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(5,6) ^ Sj s *>CRj> Uj + 0((^ u^) 2 + (Z\. Uj) 2 ), if h T (w,w> £ 0 
and/or 

= *>(Rj) C\. Uj + 0((A- u^ 2 + (^ Uj) 2 ), if h Q (w,w) £ 0 
Dropping the quadratic terms, inequality (2.9) gives 


*OQ 

(5.7) (a) 0 < — i)(RT) < 1 if h.(w,w) 4 0 

r“ J 1 


♦<0 . 

(b) 0 < -J-, l)(Q < 1 if h.(w,w> t 0 

R J 

These inequalities are a bit more restrictive than those needed in 
the monotone case. There we have 



i>aO 


+ 4KR. .) if h. (w,w) £ 0 


j m1 ) if h 0 (w,w) £ 0 

which. are the same as inequalities (5.4) (a) ,(b) . 

For the entropy condition to be proven valid, the flux was modi- 
fied. V/e replaced H by H ac (where a.c. stands for artificial 


>MRt) 

(b) 1 > V 2 “ ^ + ^( R 

Rt 

J 
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compression), as follows, 


First we modify the flux differences in (5.2): 


(5.9) (a) (h(u j+1 , u j) " f(u j )) 


* <h(u j + i- 'V - f(u j )) 1 + 


h i tu .r u j i 

h( Vr u j> - f( V 


(f( Uj+1 ) -h(u J+1 , Uj)) 


*;U,(/yVv va. u : 


' <f(U J+1 ) ' h(u J*1’ V’ I 1 ■ f<u J+t ) -h(u J+1 , Uj) 


Here a j + 1/ 2 are both positive, chosen first so that the quantities 
in the brackets in (5.9) are between 0 and 2, (inequalities (5.13) 
below) . 


We next let: 


, H + . (f(u i ) - h(u .i' Vi” 

J -«v,. V )M 


. (hl “w V -«"1» 

d (h( V Vi } - f( Vi ))M 
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Then we define our scheme: 


(5.10) 


OU . 1 

"St = - AA H (u j +2 - U j + 1 - V U J- 1 > 


(5.11) H a “(u J+2> u. +1 , « Uj .,) 


• h(u j+i> V - Tfe *«j';i> (h(u j+r V - 


+ V 2 Wf) (f(u j+1 > - h(u j+1 , Uj )) M . 


The resulting scheme is convergent if the following inequalities 


are valid: 


(5.12) X u J + 1 f«(w)dw ((V 2 (Z\ u.)) 2 - (w - V 2 (u j+1 + Uj)) 2 ) 


- Vv 2 (A, Uj) [>yUj + i. u j+1 ) - h,( Ujt Uj)] 


- a j + v 2 'A, u j > 2 lh o (u j.i>Vi> - h o ( V u j n i 0 


(h(u, 1t u.) - f(u.)) 

(5. 13) (a) ia J+V2 (h,(u j+1 . u J+ ,) - h,(u., Uj»| < 1— *- 


f(u J+1 } - h(u .1+1« u 1> 


|a W h o ( Vi’ Vi 5 "VV “A 1 j — 7^1 
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Examples are given in [12], 


It is interesting that this case allows compression to be added at 

4 * 

shocks, i.e. the restrictions on the a“ i. allow negative viscosity to 

J + >2 

be added to (5.2) if > u^, while, again for the proof of conver- 
gence, positive viscosity must be added for rarefactions, < Uj + ^. In 
the present MUSCL case the situation is (unnaturally) reversed. A com- 
parison of the convergent MUSCL scheme with this one yields two main 
points: 


(1) Away from sonic points, the fluxes differ only in 0((/^ u^)^) 


terms 


(2) At sonic points, the fluxes differ in 0((^ u^) 2 ) terms, 
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